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Abstract 

We propose a new duality involving topological strings in the limit 
of large string coupling constant. The dual is described in terms of 
a classical statistical mechanical model of crystal melting, where the 
temperature is inverse of the string coupling constant. The crystal is 
a discretization of the toric base of the Calabi-Yau with lattice length 
gs. As a strong evidence for this duality we recover the topological 
vertex in terms of the statistical mechanical probability distribution 
for crystal melting. We also propose a more general duality involving 
the dimer problem on periodic lattices and topological A-model string 
on arbitrary local toric threefolds. The {p, q) 5-brane web, dual to 
Calabi-Yau, gets identified with the transition regions of rigid dimer 
configurations. 
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1 Introduction 

Topological strings on Calabi-Yau threefolds haye been a fascinating class of 
string theories, which haye led to insights into dynamics of superstrings and 
super symmetric gauge theories. They haye also been shown to be equiyalent 
in some cases to non-critical bosonic strings. In this paper we ask how the 
topological A-model which 'counts' holomorphic curyes inside the Calabi-Yau 
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behaves in the hmit of large values of the string coupling constant ^ 1. 
We propose a dual description which is given in terms of a discrete statistical 
mechanical model of a three dimensional real crystal with boundaries, where 
the crystal is located in the toric base of the Calabi-Yau threefold, with the 
'atoms' separated by a distance of gs- Moreover the Temperature T in the 
statistical mechanical model corresponds to 1/gs- Heating up the crystal 
leads to melting of it. In the limit of large temperature, or small gs, the 
Calabi-Yau geometry emerges from the geometry of the molten crystal ! 

In the first part of this paper we focus on the simplest Calabi-Yau, namely 
C'^. Even here there are a lot of non-trivial questions to answer. In particular 
the computation of topological string amplitudes when we put D-branes in 
this background is non-trivial and leads to the notion of a topological vertex 
[H Ej (see also the recent paper [3]). Moreover using the topological vertex 
one can compute an all order amplitude for topological strings on arbitrary 
local Calabi-Yau manifolds. This is an interesting class to study, as it leads to 
non-trivial predictions for instanton corrections to gauge and gravitational 
couplings of a large class of = 2 supersymmetric gauge theories in 4 
dimensions via geometric engineering |3] (for recent progress in this direction 
see |SJin]). It is also the same class which is equivalent (in some limits) to non- 
critical bosonic string theories. In this paper we will connect the topological 
vertex to the partition function of a melting corner with fixed asymptotic 
boundary conditions. Furthermore we find an intriguing link between dimer 
statistical mechanical models and non-compact toric Calabi-Yau threefolds. 
In particular the dimer problems in 2 dimensions naturally get related to the 
study of configurations of {p, q) 5-brane web, which is dual to non-compact 
toric Calabi-Yau threefolds. 

The organization of this paper is as follows: In Section 2 we will motivate 
and state the conjecture. In Section 3 we check aspects of this conjecture 
and derive the topological vertex from the statistical mechanical model. In 
Section 4 we discuss dimer problems and its relation to topological strings 
on Calabi-Yau. 

Our proposal immediately raises many questions, which are being presently 
investigated. Many of them will be pointed out in the paper. One of the 
most interesting physical questions involves the superstring intepretation of 
the discretization of space. On the mathematical side, we expect that our 
statistical mechanical model should have a deep meaning in terms of the 
geometry of target space based on the interpretation of its configurations as 
torus fixed points in the Hilbert scheme of curves of the target threefold. 
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Also, our 3d model naturally extends the random 2d partition models that 
arise in = 2 supersymmetric gauge theory j3] and Gromov-Witten theory 
of target curves [7] . For some mathematical aspects of the topological vertex 
see IHliniE!. 



2 The Conjecture 

2.1 Hodge integrals and 3d partitions 

Consider topological A-model strings on a Calabi-Yau threefold. For simplic- 
ity, let us consider the limit when the Kahler class of Calabi-Yau is rescaled 
by a factor that goes to infinity. As explained in jTHj in this limit the genus 
g amplitude is given by 



X 



where x is the Euler characteristic of the Calabi-Yau threefold, is the 
moduli space of Riemann surfaces, 5{ is the Hodge bundle over Mg, and Cg-i 
denotes the {g — l)-st Chern class of it. For genus and 1, there are also 
some Kahler dependence (involving volume and the second Chern class of 
the tangent bundle), which we subtract out to get a finite answer. Consider 



Z = exp 



Eds'-' i 4-m 



It has been argued physically [TT] and derived mathematically [12] that 

Z = 

where 

/=n t 

■> ii(l_gn)n 

n ^ ' 

and 

q = e"^" . 

By the classical result of McMahon, the function / is the generating function 
for 3d partitions, that is. 



# boxes 

3d partitions 



where, by definition, a 3d partition is a 3d generalization of 2d Young dia- 
grams and is an object of kind seen on the left in Figure 01 This fact was 
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pointed out to one of US cLS cL curiosity by R. Dijkgraaf shortly after [TT] 
appeared. 

2.2 Melting of a crystal and Calabi-Yau threefold 

It is natural to ask whether there is a deeper reason for this correspondence. 
What could three dimensional partitions have to do with A-model topological 
string on a Calabi-Yau threefold ? The hint comes from the fact that we 
are considering the limit of large Kahler class and in this limit the Calabi- 
Yau looks locally made of C^'s glued together. It is then natural to view 
this torically, as we often do in topological string, in the context of mirror 
symmetry, and write the Kahler form as 



and \zi\^ span the base of a toric fibration of C^. Note that \zi\^ = Xi 
parameterize the positive octant O"*" C R'^. If we assign Euler characteristic 
"2" to each patch, the topological string amplitudes on it get related to 
the McMahon function. Then it is natural to think that the octant is related 
to the three dimensional partitions, in which the boxes are located at 7? 
lattice points inside 0"*~. Somehow the points of the Calabi-Yau, in this case 
C^, become related to integral lattice points on the toric base. 

The picture we propose is the following. We identify the highly quantum 
Calabi-Yau with the frozen crystal, that is, the crystal in which all atoms 
(indexed by lattice points in O^) are in place. We view the excitations as 
removing lattice points as in Figure H The rule is that we can remove lattice 
points only if there are no pairs of atoms on opposite sides. This gives the 
same rule as 3d partitions. Note that if one holds the page upside-down, one 
sees a 3d partition in Figure ^ namely the partition from Figure El 

Removing each atom contributes the factor q = e~^^'^ to the Boltzmann 
weight of the configuration, where fi is the chemical potential (the energy 
of the removal of an atom) and T is the temperature. We choose units in 
which /i = 1. To connect this model to topological string on we identify 
Qs = 1/T. In particular, the Qs — >■ 0, that is, the g — >■ 1 limit the crystal 
begins to melt away. Rescaled in all direction by a factor of 1/T = Qs, the 
crystal approaches a smooth limit shape which has been studied from various 
viewpoints ^] . It is plotted in Figure El 

The analytic form of this limit shape is encoded in terms of a complex 




i=l,2,3 



j=l,2,3 
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Figure 1: A melting crystal corner 




5 5 

Figure 2: The limit shape of a 3d partition 
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Riemann surface, which in this case is given by 



t;) = e"" + e"'' + 1 = (2.1) 

defined as a hypersurface in with a natural 2-form du A dv, where u, v 
are periodic variables with period 27rz. In coordinates e~" and this is 
simply a straight line in C^. Consider the following function of the variables 
U = Re{u) and V = Re{v) 

R{U,V) = ^jj^ log \F{U + i9,V + i4>)\d9d^. (2.2) 

This function is known as the Ronkin function of F. In terms of the Ronkin 
function, the limit shape can be parameterized as follows: 

{xi,X2,X3) = {U + R,V + R,R), R = R{U,V). 

Note that 

U = Xi - X3 , V = X2 - X-i . 

The projection of the curved part of the limit shape onto the (f/, V^)-plane is 
the region bounded by the curves 

±e-^ ± e"^ + 1 = , 

excluding the case when both signs are positive. This region is the amoeba 
of the curve ()2.1|) . which, by definition, is its image under the map {u,v) ^ 
{U,V). In different coordinates, this is the planar region actually seen in 
Figure |21 



2.3 Mirror symmetry and the limit shape 

Consider topological strings on C^. One can apply mirror symmetry in this 
context by dualizing the three phases of the complex parameters according 
to T-duality. This has been done in [12] . One introduces dual variables Yi 
which are periodic, with period 27ri and are related to the Zi by 

\zi\^ = ReYi 

(this is quantum mechanically modified by addition of an ^-independent large 
positive constant to the right hand side). The imaginary part of the Yi is 
'invisible' to the original geometry, just as is the phase of the Zt to the Y^ 
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variables. They are T-dual circles. But we can compare the data between 
and the mirror on the base of the toric variety which is visible to both. 
In mirror symmetry the Kahler form of C'^ gets mapped to the holomorphic 
three form which in this case is given by 

3 
i=l 

where 



Note that shifting — > Kj + r shifts W — > e~^W. The analog of rescaling 
in the mirror is changing the scale of W. Let us fix the scale by requiring 
W = 1; this will turn out to be the mirror statement to rescaling by Qs to 
get a limit shape. In fact we will see below that the limit shape corresponds 
to the toric projection of the complex surface in the mirror given simply by 
WiYi) = 1. Let us define 

u = Y^-Y,, v = Y,-Y3 

Then 

W = e-^'F{u, v) 

where 

F{u, v) = e"" + e"'' + 1 

we will identify F{u,v) with the Riemann surface of of the crystal melting 
problem. In this context also according to mirror map the points in the u, v 
space get mapped to the points on the toric base (i.e 0~^) satisfying 

xi — X3 = Re{u) = U 

X2 — X3 = Re{v) = V 

We now wish to understand the interpretation of limit shape from the view- 
point of topological string. We propose that the boundary of the molten 
crystal which is a 2-cycle on the octant should be viewed as a special La- 
grangian cycle of the A-model with one hidden circle in the fiber. Similarly in 
the B-model mirror it should be viewed as the B-model holomorphic surface, 
which in the case at hand gets identified with W = 1 (recall that in the LG 
models B-branes can be identified with W = const. ^1]). On this surface we 
have 

e-^'F{u,v) = 1 
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If we take the absolute value of this equation, to find the projection on to 
the base we find 

If we take the logarithm of this relation we have 

-ReiYs) + log|e"" + + 1| = ^ 0:3 = log|e~" + e^'' + 1| 

However we have a fuzziness in mapping this to the toric base: u = Re{u) +i6 
and V = Re{v) + i(f) and so a given value for Re{u) and Re{v) does not give a 
fixed value of X3. That depends in addition on the angles 6*, (f) of the mirror 
torus which are invisible to the A-model toric base. It is natural to take the 
average values as defining the projection to the base, i.e. 

X3 = ^ y dOdct) \og\F{U + ie,V + i(t))\ 

which is exactly the expression for the limit shape. We thus find some further 
evidence that the statistical mechanical problem of crystal melting is rather 
deeply related to topological string and mirror symmetry on Calabi-Yau. 
Moreover we can identify the points of the crystal, with the discretization of 
points of the base of the toric Calabi-Yau. 

To test this conjecture further we will have to first broaden the dictionary 
between the two sides. In particular we ask what is the interpretation of the 
topological vertex for the statistical mechanical problem of crystal melting? 
For topological vertex we fix a 2d partition on each of the three legs of the 
toric base. There is only one natural interpretation of what this could mean 
in the crystal melting problem: This could be the partition function of the 
melting crystal with three fixed asymptotic boundary shapes for the molten 
crystal, dictated by the corresponding partition. We will show this is indeed 
the case in the next section. 



3 Melting corner and the topological vertex 
3.1 Transfer matrix approach 

The grand canonical ensemble of 3d partitions weighted with g* is the 
simplest model of a melting crystal near its corner. We review the transfer 
matrix approach to this model following ^Hj- This approach can be easily 
generalized to allow for certain inhomogeneity and periodicity, which is useful 
in the context of more general models discussed in Section IH 
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We start by cutting the 3d partition into diagonal slices by planes X2—X1 = 
t, see Figure El This operation makes a 3d partition a sequence {/i(t)} of 
ordinary partitions indexed by an integer variable t. Conversely, given a 
sequence it can be assembled in a 3d partition provided it satisfies 

the following interlacing condition. We say that two partitions fi and u 
interlace, and write fi >~ u ii 

fJ'i > J^i > 1^2 > 1^2 > ■ ■ ■ ■ 
It is easy to see that a sequence of slices of a 3d partition satisfies 

fi{t) ^ fi{t + 1) , t<0, 
and the reverse relation for t > 0. 




Figure 3: A 3d partition and its diagonal slices 



There is a well-known map from partitions to states in the NS sector of 
complex fermionic oscillator. Let and bi be the areas of pieces one gets by 
slicing a 2d partition first diagonally and then horizontally (resp. vertically) 
above and the below the diagonal, respectively. Formally, 

ai = fj.i-i + -, bi = + 

where i ranges from 1 to the number of squares on the diagonal of fi. In math- 
ematical literature, these coordinates on partitions are known as (modified) 
Frobenius coordinates. The fermionic state associated to fi is 

d 
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Note that 

qLo\^) =q# =gl^'|/i), 

where we denote the total number of boxes of /i by We can write a 
bosonic representation of this state by the standard bosonization procedure. 

Consider the operators 



±n>0 / 



where J„ denotes the modes of the fermionic current The operators 

T±{z) can be identified with annihilation and creation parts of the bosonic 
vertex operator e'^^^\ The relevance of these operators for our problem lies 
in the following formulas: 

r-(i)l/i) = El^) (3.3) 
r+(i)|/x) = 5^k) 



To illustrate their power we will now derive, as in ^H]; the McMahon's 
generation function for 3d partitions 

^ ^ ^ g# of boxes 

3d partitions tt 

By the identification ()3.3p of the transfer matrix, we have 

^=((n^'°r+(l)jg^"(^n T^{l)q^^Y (3.4) 

Now we commute the operators to the outside, splitting the middle one 
in half. This yields 

^=(nr+(5"^^)nr-(9^""^))- (3.5) 

\n.>0 n>0 / 

Now we commute the creation operators through annihilation operators, us- 
ing the commutation relation 

T^{z)T.{z') = {l-z/zr'T^{z')T^{z). 
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The product of resulting factors gives directly the McMahon's function 

z = i[{i-q-r^ = f. 

n>0 

Same ideas can be used to get more refined results, such as, for example, 
the correlation functions, see ^H] • We will generalize below the above compu- 
tation to the case when 3d partitions have certain asymptotic configuration 
in the direction of the three axes. Before doing this, we will review some 
relevant theory of symmetric functions. 

3.2 Skew Schur functions 

Skew Schur functions Sx/^{xi, X2, ■ ■ ■) are certain symmetric polynomials in 
the variables Xi indexed by a pair of partitions /i and A such that /i C A, see 
P!7] . Their relevance for us lies in the well-known fact (see e.g. ^H]) that 

llT.ix,)\fi) = J2^x/,i^)\>^)- (3-6) 

i XD/j, 

When /i = 0, this specializes to the usual Schur functions. The Jacobi- Trudy 
determinantal formula continues to hold for skew Schur functions: 

sa/m = det {hx,-t,j+j-i) . (3.7) 

Here hk is the complete homogeneous function of degree k — sum of all 
monomials of degree k. They can be defined by the generating series 

J2hJ^ = l[{l-tx,r'=exp(Y.^-^Y.^q ■ (3.8) 

n>0 i \n>0 i / 

The formula (j3.7|) is very efficient for computing the values of skew Schur 
function, including their values at the points of the form 

= (g-i-V2^ ^.2-3/2^ ^.3-5/2^ _ _ ^ ^ (3 

where z/ is a partition. In this case the sum over i in ()3.8j) becomes a geometric 
series and can be summed explicitly. 

There is a standard involution in the algebra of symmetric function which 

acts by 

s\ ^ Sxt , 
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where A* denotes the transposed diagram. It continues to act on skew Schur 
functions in the same manner 



It is straightforward to check that 

sx/,{q''^') = i-iy'^-^'^s^y^.iq-'-^) . (3.10) 

Finally, the following property of the skew Schur function will be crucial 
in making the connection to the formula for the topological vertex from 
The coefficients of the expansion 



sx/^, = ^4uSu (3.11) 



1/ 



of skew Schur functions in terms of ordinary Schur function are precisely 
the tensor product multiplicities, also known as the Littlewood-Richardson 
coefficients. 

3.3 Topological vertex and 3d partitions 
3.3.1 Topological vertex in terms of Schur functions 

Our goal now is to recast the topological vertex P in terms of Schur functions. 
The basic ingredient of the topological vertex involves the expectation values 
of U{oo) Chern-Simons Hopf link invariant VF^a = Wx^ in representations fi 
and A. The expression of the Hopf link invariant W^^x in terms of the Schur 
functions is the following: 



where 
Here 



W,x = W^sx{q^+n, (3.12) 
W^;. = g''^^)/'V(/)- (3-13) 



K 

' / J I \ " 'A I \ Z / 

□=(i,j)GA 

is the unique up to scalar quadratic Casimir such that 

ft(0) = ft(n) =0. 
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Using P-lip it is straightforward to check that the topological vertex C(A, fi, v) 
in the standard framing has the following expression in terms of the skew 
Schur functions 

C(A,/i,z/) =g''W/2+'^('^)/2s..(gO s.v.lg'^^O • (3-15) 

3.3.2 The lattice length 

To relate the crystal melting problem to the topological vertex we first have 
to note that the topological vertex refers to computations in the A-model 
corresponding to placing Lagrangian branes on each leg of C^, assembled 
into representations of f/(oo) and identified with partitions. If we place the 
brane at a fixed position and put it in a representation /x, the effect of moving 
the brane from a position / to the position I + k (in string units) affects the 
amplitude by a multiplication of 

exp(— /c|/i|) 

Now consider the lattice model where we fix the asymptotics at a distance 
L >> 1 to be fixed to be a fixed partition /x. Then if we change L ^ L + K 
then the amplitude gets weighted by g^l^L Since we have already identified 
q = e"^" Comparing these two expressions we immediately deduce that 

k = Kqs 

In other words the distance in the lattice computation times Qs is the distance 
as measured in string units. This is satisfactory as it suggests that as Qs —>■ 
the lattice spacing in string units goes to zero and the space becomes 
continuous. 

In defining the topological vertex one gets rid of the propagator factors 
above (which will show up in the gluing rules). Similarly in the lattice model 
when we fix the asymptotic boundary condition to be given by fixed 2d par- 
titions we should multiply the amplitudes by q~^^f^^ for each fixed asymptote 
at lattice position L. Actually this is not precisely right: we should rather 
counter weight it with ^"(-^+2)1^1. To see this note that if we glue two topolog- 
ical vertices with lattice points Li and L2 along the joining edge, the number 
of points along the glued edge is Li + L2 + 1. Putting the | in the above 
formula gives a symmetric treatment of this issue in the context of gluing. 
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3.3.3 Framing 

Topological vertex also comes equipped with a framing P for each edge, 
which we now recall. Toric Calabi-Yau's come with a canonical direction in 
the toric base. In the case of C"^, it is the diagonal line Xi = X2 = X3 in the 
O^. One typically projects vectors on this toric base along this direction, 
to a 2-dimensional plane (the U, V plane in the context we have discussed). 
At each vertex there are integral projected 2d vectors along the axes which 
sum up to zero. The topological vertex framing is equivalent to picking a 2d 
projected vector on each axis whose cross product with the integral vector 
along the axis is +1 (with a suitable sense of orientation). If Vi is a framing 
vector for the i-th axis, and Cj denotes the integral vector along the i-th axis, 
then the most general framing is obtained by 

where rii is an integer. We now interpret this choice in our statistical me- 
chanical model: In describing the asymptotes of the 3d partition we have to 
choose a slicing along each axis. We use the framing vector, together with 
the diagonal direction xi = X2 = X3 to define a slicing 2-plane for that edge. 
The standard framing corresponds to choosing the framing vector in cyclic 
order: On the xi-axis, we choose 0:3, on the X3-axis we choose X2 and on the 
X2-axis we choose Xi. This together with the diagonal line determines a slic- 
ing plane on each axis. Note that all different slicings will have the diagonal 
line on them. This line passes through the diagonal of the corresponding 2d 
partition. 

Before doing any detailed comparison with the statistical mechanical model 
with fixed asymptotes we can check whether framing dependence of the topo- 
logical vertex can be understood. This is indeed the case. Suppose we 
compute the partition function with fixed 2d asymptotes and with a given 
framing (i.e. slicing). Suppose we shift the framing by n^. This will still 
cut the asymptotic diagram along the same 2d partition. Now, however, the 
total number of boxes of the 3d partition has changed. The diagonal points 
of the partition have not moved as they are on the slicing plane for each 
framing. The farther a point is from the diagonal the more it has moved. 
Indeed the net number of points added to the 3d partition is given by 

ni'^ik - I) = riiK^fi) 

Thus the statistical mechanical model will have the extra Boltzmann weight 
qUiKiiJ.) xhis is precisely the framing dependence of the vertex. Encouraged 
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by this observation we now turn to computing the topological vertex in the 
standard framing from the crystal point of view. 

3.4 The perpendicular partition function 
3.4.1 Definition 

Now our goal is to find an exact match between the formula ()3.15|) and 
the partition function P(A, /x, u) for 3d partitions whose asymptotics in the 
direction of the three coordinates axes is given by three given partitions A, 
fi, and u. This generating function, which we call the perpendicular partition 
function, will be defined and computed presently. 
Consider 3-dimensional partitions vr inside the box 

[0,N^] X [CiVj X [0,N,]. 

Let the boundary conditions in the planes Xi = Ni be given by three partitions 
A, /i, and v. This means that, for example, the facet of vr in the plane Xi = Ni 
is the diagram of the partition A oriented so that Ai is its length in the X2 
direction. The other two boundary partitions are defined in the cyclically 
symmetric way. See FigureEJin which A = (3, 2), /i = (3, 1) and u = (3, 1, 1). 



Figure 4: A 3d partition ending on three given 2d partitions 



Let Pni,N2,N3 (a, yU, i^) be the partition function in which every n is weighted 
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by It is obvious that the hmit 

P(A,/i,i.) = hm g-^^l^l-^^H-^3klp^^^^^^^^(A,^,^) (3.I6) 

^1,^2,^3— >00 

exists as a formal power series in q. What should the relation of this to 
topological vertex be? From what we have said before this should be the 
topological vertex itself up to framing factors. Let us also fix the framing 
factor, to compare it to the canonical framing. First of all we need to multiply 
P by 

This is related to our discussion of the gluing algorithm (of shifting Ni Ni+ 
^) and splitting the point of gluing between the two vertices. Secondly this 
perpendicular slicing is not the same as canonical framing. We need to rotate 
the perpendicular slicing to become the canonical framing. This involves the 
rotation of the partition along its first column by one unit and the sense of 
the rotation is to increase the number of points. For each representation this 
gives 

^5^A.(A.-l) = ^(||Af -|A|) 

i 

extra boxes which we have to subtract off and gives the additional weight. 
Here 

Combining with the previous factor we get a net factor of 

g|(l|A|P + ||MlP+llHl') 

Moreover we should normalize as usual by dividing by the partition function 
with trivial asymptotic partition, which is the McMahon function. We thus 
expect 

- g")"gi(ll^lP + IHP+IIHP)p(A, ^, u) = C(A, 12, V). 

n 

We will see below that this is true up to Qs —Qs and an overall factor that 
does not affect the gluing properties of the vertex (as they come in pairs) 
and it can be viewed as a gauge choice for the topological vertex. 

3.4.2 Transfer matrix formula 

Recall that in the transfer matrix setup one slices the partition diagonally. 
Compared with the perpendicular cutting, the diagonal cutting adds extra 
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boxes to the partition and, as a result, it increases its volume by Q) + ('2)' 
where, by definition 

Also, for the transfer-matrix method it is convenient to let N3 = 00 from the 
very beginning, that is, to consider 

PNuN2{\fJ',^) = lim q~^''^''^PNuN2,N3{\ ■ (3.18) 

A^s— >oo 

The partition function Pn^^n^IX, fi, v) counts 3d partitions vr with given bound- 
ary conditions on the planes Xi^2 = Ni^2 inside the container which is a semi- 
infinite cylinder with base 

[0,iVi] X [0,iV2]\i/, 

see Figure EJ In other words, Pni,N2{X, fi, u) counts skew 3d partitions in the 
sense of |TBj . 




Figure 5: The container for skew 3d partitions 



The main observation about skew 3d partitions is that their diagonal slices 
interlace in the pattern dictated by the shape v. This gives the following 
transfer-matrix formula for ()3.18jl . which is a direct generalization of ()3.4|1 : 



\ [ 



I terms 



r±(i) 



n 

N2-1 

^terms 



r±(i) 



fi ) , (3.19) 
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where the pattern of pluses and minuses in T± is dictated by the shape of u. 

3.4.3 Transformation of the operator formula 

Now we apply to the formula ()3.19p the following three transformations: 

• Commute the operators to the outside, splitting the middle one in 
half. The operators r-|-(l) will be conjugated to the operators r-|-(g---) 
in the process. 

• Commute the raising operators r_ to the left and the lowering operators 
r+ to the right. There will be some overall, i^-dependent multiplicative 
factor Z^u) from this operation. 

• Write the resulting expression as a sum over intermediate states \ri){ri\. 

3.4.4 The result 

We obtain 



P(A, u) = Z{u) q~{l)<iy\'\l^~\^\l^ Yl ^AV.(?-^-0 s,/,{q-'^'-n ■ (3.20) 



V 



In order to determine the multiplicative factor Z{u), we compute 



P(0, 0, i^) = P(z/, 0, 0) 



using the formula ()3.20|) . We get 



-(^)H-l/^..(g-^) 

nn>o(i-^")'^ 



(3.21) 



Since 




(3.22) 



comparing (|3.2(Jj) with (j3.15|) we obtain 



C(A,/i, z/; 1/q) = q 



ll>-*|| + ||M'll + lk*ll 



l[{l-qn''P{X,fi,u) 



(3.23) 



2 



n>0 



where 




(3.24) 
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This is exactly what we aimed for, see Section 13.4.11 up to the following 
immaterial details. The difference between ||A*|| and ||A|| is irrelevant since 
the gluing formula pairs A with A*. Also, since the string amplitudes are 
even functions of (?s, they are unaffected by the substitution g i— s> 1/g. These 
differences can be viewed as a gauge choice for the topological vertex. 

3.5 Other generalizations 

One can use the topological vertex to glue various local patches and obtain 
the topological A-model amplitudes. Thus it is natural to expect that there 
is a natural lattice model. There is a natural lattice [22] for the crystal in this 
case, obtained by viewing the Kahler form divided by Qs as defining the first 
Chern class of a line bundle and identifying the lattice model with holomor- 
phic sections of this bundle. This is nothing but geometric quantization of 
Calabi-Yau with the Kahler form playing the role of symplectic structure and 
Qs playing the role of h. The precise definition of lattice melting problem for 
this class is currently under investigation [2^1 • The fact that we have already 
seen the emergence of topological vertex in the lattice computation corre- 
sponding to one would expect that asymptotic gluings suitably defined 
should give the gluing rules of the statistical mechanical model. It is natural 
to expect this idea also works the same way in the compact case, by view- 
ing the Calabi-Yau as a non-commutative manifold with non-commutativity 
parameter being Qsk, where k is the Kahler form. 

It is also natural to embed this in superstring (this is currently under 
investigation j2ZI), where gs will be replaced by graviphoton field strength. 
The large in this context translates to strong graviphoton field strength, 
for which it is natural to expect discretization of spacetime. This is exciting 
as it will potentially give a novel realization of superstring target space as a 
discrete lattice. 

4 Periodic dimers and toric local CY 

4.1 Dimers on a periodic planar bipartite graph 

A natural generalization of the ideas discussed here is the planar dimer model, 
see fH] for an introduction. Let F be a planar graph ("lattice") which is 
periodic and bipartite. The first condition means that it is lifted from a 
finite graph in the torus via the standard covering map T^. The 
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second condition means that the vertices of T can be partitioned into two 
disjoint subsets ("black" and "white" vertices) such that edges connect only 
white vertices to black vertices. Examples of such graphs are the standard 
square or honeycomb lattices. 

By definition, a dimer configuration on F is a collection of edges D = {e} 
such that every vertex is incident to exactly one edge in D. Subject to 
suitable boundary conditions, the partition function of the dimer model is 
defined by 

D eeD 

where w{e) is a certain (Boltzmann) weight assigned to a given edge. For 
example, one can take both weights and boundary conditions to be periodic, 
in which case Z is a finite sum. One can also impose boundary conditions at 
infinity by saying that the dimer configuration should coincide with a given 
configuration outside some ball of a large radius. The relative weight of such 
a configuration is a well-defined finite product, but the sum Z itself is infinite. 
In order to make it convergent, one introduces a factor of g™iumc(D)^ defined 
in terms of the height function, see below. Other boundary conditions can be 
given by cutting a large but finite piece out of the graph F and considering 
dimer s on it. 

Simple dimer models, such as equal weight square or honeycomb grid dimer 
models with simple boundary conditions were first considered in the physics 
literature many years ago (201 A complete theory of the dimer model 
on a periodic weighted planar bipartite graph was developed in |211I2S1- It 
has some distinctive new features due to spectral curve being a general high 
genus algebraic curve. We will now quote some results of [2111221 and indicate 
their relevance in our setting. 

4.2 Periodic configurations and spectral curve 

Consider dimer configuration D on a torus or, equivalently, a dimer configu- 
ration in the plane that repeat periodically, like a wall-paper pattern. There 
are finitely many such configurations and they will play a special role for us, 
namely they will describe the possible facets of our CY crystal. We will now 
introduce a certain refined counting of these configurations. 

Given two configurations Di and D2 on a torus T^, their union is a collec- 
tion of closed loops on T^. These loops come with a natural orientation by, 
for example, going from white to black vertices along the edges of the first 
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dimer and from black to white vertices along the edges of the second dimer. 
Hence, they define an element (by summing over all classes of the loops) 

h = {h,h2)eH^{T\Z), 

of the first homology group of the torus. Fixing any configuration Dq as our 
reference point, we can associate h = h{D) to any other dimer configuration 
and define 

F{z, w) = ^(-l)QW w''^ W w{e) , (4.25) 

D eGD 

where Q{h) is any of the 4 theta-characteristic, for example, Q{{hi, h2)) = 
hih2- The ambiguity in the definition of ()4.25|) comes from the choice of the 
reference dimer Dq, which means overall multiplication by a monomial in z 
and w, the choice of the basis for iJi(T^, Z), which means SL{2, Z) action, 
and the choice of the theta-characteristic, which means flipping the signs of 
z and w. 

The locus F{z, w) = defines a curve in the toric surface corresponding to 
the Newton polygon of F. It is called the spectral curve of the dimer problem 
for the given set of weights. It is the spectral curve of the Kasteleyn operator 
on r, the variables z and w being the Bloch-Floquet multipliers in the two 
directions. In our situation, the curve F{z,w) will be related to the mirror 
Calabi-Yau threefold. 

4.3 Height function and "empty" configurations 

Now consider dimer configurations in the plane. The union of two dimer 
configurations Di and D2 again defines a collection of closed oriented loops. 
We can view it as the boundary of the level sets of a function h, defined on 
the cells (a.k.a. faces) of the graph F. This function h is well-defined up to 
a constant and is known as the height function. 

It is instructive to see how for dimers on the hexagonal lattice this re- 
produces the combinatorics of the 3d partitions, the height function giving 
the previously missing 3rd spatial coordinate, see Figure IHl The "full cor- 
ner" or "empty room" configuration, which was our starting configuration 
describing the fully quantum in the language of the dimers becomes the 
unique, up-to translation, configuration in which the periodic dimer patterns 
can come together. Each rhombus corresponds to one dimer (the edge of the 
honeycomb lattice inside it). 
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Figure 6: The "Empty room" configuration of honeycomb dimers 



For general dimers, there are many periodic dimer patterns and there 
are (integer) moduli in how they can come together. For example, for the 
square lattice, which is the case corresponding to the 0(— 1) © 0(— 1) — ^ 
geometry, the periodic patterns are the brickwall patterns and there is one 
integer degree of freedom in how they can be patched together. One possible 
such configuration is shown in Figure [7| The arrows in Figure [7| point from 
white vertices to black ones to help visualize the difference between the 4 
periodic patterns. 

All of these "empty" configurations can serve as the initial configuration, 
describing the fully quantum toric threefold, for the dimer problem. When 
lifted in 3d via the height function, each empty configuration follows a piece- 
wise linear function, which is the boundary of the polyhedron defining the 
toric variety. In particular, the number of "empty room" moduli matches 
the Kahler moduli of the toric threefold and changes in its combinatorics 
correspond to the flops. 

Readers familiar with [p, q) 5-brane webs [221 and their relation to toric 
Calabi-Yau j2S] immediately see a dictionary: A 5-brane configuration is 
identified as the transition line from one rigid dimer configuration to other, 
which in toric language is related to which cycle of Calabi-Yau shrinks over 
it. For example the the geometry is mapped to three 5-branes (of types 
(1, 0), (0, 1), (—1, —1)) on the 2 plane meeting at a point where each region 
gets identified with a particular dimer configuration. Similar description 
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Figure 7: An "Empty room" configuration of square dimers 



fiolds for arbitrary 5-brane webs. In tliis context the F{z,w) = F{e ",e ^) 
is identified witli tlie mirror geometry P^.^ 

4.4 Excitations and limit shape 

Now we can start "removing atoms" from tlie "full crystal" configuration, 
adding the cost of q to each increase in the height function. The limit shape 
that develops, is controlled by the surface tension of the dimer problem. This 
is a function of a slope measuring how much the dimer likes to have height 
function with this slope. Formally, it is defined as the n ^ oo limit of the 
free energy per fundamental domain for dimer configurations on the n x n 
torus restricted to lie in a given homology class or, equivalently, restricted 
to have certain slope when lifted to a periodic configuration in the plane. 

One of the main results of is the identification of this function with 
the Legendre dual of the Ronkin function of the polynomial ()4.25p defined 



The Wulff construction implies that the Ronkin function itself is one of the 
possible limit shapes, the one corresponding to its own boundary conditions. 
Note that this is exactly what one would anticipate from our general con- 
jecture if we view F{z, w) as the describing the mirror geometry. In fact 
following the same type of argument as in the C"^ case discussed before, 
would lead us to the above Ronkin function. 

As an example consider the dimers on the hexagonal lattice with 1x1 
fundamental domain. In this case the edge weights can be gauged away and 

^More precisely as in |13lll4j this corresponds to a LG theory with W = e^^'(F(u,w)) 
or to a non-compact CY given as a hypersurface: a/? — F{u, v) = 0. 



by 




•2: 



(4.26) 
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()4.26|) becomes the Ronkin function considered in Section For the square 
lattice with 1x1 fundamental domain, there is 1 gauge invariant combination 
of the 4 weights and the spectral curve takes the form 



where t is a parameter related to the size of in the 0(— 1)©0(— 1) — s>P^ 
geometry. The spectral curve is a hyperbola in and (the negative of) its 
Ronkin function is plotted in Figure |H1 Note how one can actually see the 
projection of the mirror curve ! 



All possible limit shapes for given set of dimer weights are maximizers of 
the surface tension functional and in this sense are very similar to minimal 
surfaces. An analog of the WeierstraB parameterization for them in terms of 
analytic data was found in jSSj- It reduces the solution of the Euler-Lagrange 
PDE's to solving equations for finitely many parameters, essentially finding a 
plane curve of given degree and genus satisfying certain tangency and periods 
conditions. 

In the limit of extreme weights and large initial configurations, the amoebas 
and Ronkin functions degenerate to the piecewise-linear toric geometry. In 
this limit, it is possible to adjust parameters to reproduce the topological 
vertex formula for the GW invariants of the toric target obtained in p. We 
expect that the general case will reproduce the features of the background 
dependence (i.e. holomorphic anomaly) in the A-model fH]. This issue is 
presently under investigation fIE\ . 



F{z, w) = l + z + w — e ^zw , 




Figure 8: The Ronkin function of a hyperbola 
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